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Project Summary 


Purpose of Research 

The purpose of this research was to develop a prototype liquid hydrogen boil-off 
recovery system. Perform analyses to finalize recovery system cycle, design 
detail components, fabricate hardware, and conduct sub-component, 
component, and system level tests leading to the delivery of a prototype system 

Research Conducted and Findings 


Design point and off-design analyses, along with a sensitivity study, identified 
cycle improvements to increase the robustness of the Phase I system by adding 
a by-pass heat exchanger All components, including a high-speed, miniature 
turbocompressor, were designed and manufacturing drawings were created. All 
hardware was fabricated and tests were conducted in air, helium, and hydrogen. 
Testing validated the design and analyses, except for the turbocompressor A 
rotor-to-stator clearance issue was discovered as a result of a concentricity 
tolerance stack-up. 

Potential Applications 

A ground-based cryogenic boil-off recovery system coula be used on NASA 
propellant storage tanks. A flight-weight system could be used to minimize 
propellant losses for current launch vehicles, the space station, interplanetary 
missions, and long duration satellites. Non-government applications include 
cryogenic storage tanks at industrial gas factories, distribution centers and 
hospitals A “flight-weight” version could be used on roadable dewars used to 
transport cryogens throughout the country. 
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Project Objectives 


During Phase II, our primary objective was to develop a prototype liquid 
hydrogen boil-off recovery system and place it into service on a Stennis Space 
Center liquid hydrogen storage tank for demonstration and evaluation. This 
overall objective has been sub-divided into the following specific objectives: 

• Detail Part Design 

• Final Analysis 

• Hardware Procurement 

• Component Tests 

• System Tests 

• Prototype Tests 

Detail Design 


Create manufacturing drawings and identify off-the-shelf (OTS) components for 
all hardware that comprises the liquid hydrogen boil-off recovery system. The 
major objectives for this effort were the detail design of the compressor and 
turbine, including the blade geometry, as well as conducting the mechanical 
design necessary to generate the drawings for the unique hardware. 

Final Analysis 

Using the cycle and design point performance analysis conducted in Phase I, 
the objective of this effort was to conduct off-design and transient analysis to 
ensure the recovery system will function over a wide range of operating 
conditions, in addition to the laboratory design point. Additionally, this analysis 
allows the operating envelope, capabilities, and limitations of the recovery 
system to be identified. 

Hardware Procurement 


The focus of this objective is to acquire all hardware necessary to conduct sub- 
component, component, system, and prototype testing. Additionally, this 
objective will be used to identify manufacturing issues and refine low-volume 
production cost estimates. 

Component Tests 

This objective is to test each part at the lowest possible level beginning with sub- 
components and continuing with components. This approach will identify design 
or performance issues at the lowest possible level, which improves the 
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probability of success at higher level system and prototype tests. The bearings, 
regulator, and Jouie-Thompson valves will be tested at the sub-component level 
The turbo-compressor, heat exchangers, and separator will be tested at the 
component level. 

System Tests 

The components will be assembled into a recovery system. Testing of this 
system will be conducted with the objective of measuring operating points 
throughout the cycle. Data will be used to quantify and optimize performance at 
design point. Off-design point testing will be conducted to quantify the operating 
envelope. 

Prototype Tests 

The results from the system tests will be used to identify the prototype 
configuration and performance. This configuration will be built, tested, and 
delivered to Stennis Space Center for continued testing and field evaluation in 
accordance with the primary objective. 
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Work Performed 


Detail Design 


The detail design effort consisted of identifying all the hardware that comprises 
the liquid hydrogen boil-off recovery system. If a suitable off-the-shelf 
component could not be found, the part was designed and a manufacturing 
drawing was created. Table 1 lists the resulting bill-of-material for the recovery 
system. 

There are 37 unique parts (part numbers, not total number of pieces) that make- 
up the recovery system (note that AT9801 1 and AT98012 are assembly 
drawings and not parts and AT98010 is an alternate bearing configuration). Of 
these 38 parts (including the alternate bearing configuration), 30 are off-the- 
shelf components and 8 were custom designed by AFAB Technologies (A-Size 
versions of the detail assembly and manufacturing drawings are included as 
Appendix A). 

The eight custom parts are AT98003 - AT98010. AT98009 is the gas/liquid 
separator used in the recovery system, the remaining custom parts are used in 
the turbocompressor. Level 1 represents the entire liquid hydrogen boil-off 
recovery system. Level 2 are all th'^ oarts that make-up the recovery system 
Level 3 are the parts that make-up the turbocompressor assembly, and Level 4 
are the parts that make-up the bearing assembly. 

The majority of the design effort was spent on the turbocompressor. As such, 
the major design features of each custom part will be discussed. 
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Level 

Part Number 

1 

AT98012 1 

2 

AT98011 ( 

3 

AT98003 f 

3 

AT98004 i 


Part Name 


Description 



D12-2 


SC2-51 


SPR9-1 


AT98005 


AT98006 


AT98007 


AT98008 


AT98010 


6757782 


1384-129-01-01 


48002-3-1 5-2-SPA 


48005-5-1 5-2-SPA 


|MS9193:03 


MS91 93-05 


|A/R 


1413 


4ALOK-316 


4ET4-316 


4FSC4N-316 


4MSC4N-316 


4MSC8N-316 


4NU4-316 


I4SC4-316 


A/R 


A32 


CRMS0063035A 


ECR4-316 


SS-400-R-3 


SS-400-R-5 


A/R 


|A/R 


A/R 





Housing, Main 


Cover, Turbo com pressor 


Ring, Turbine Nozzle 


Bearing, Thrust 


Bearing, Tilting Pad 

Bolt 


Face Seal, Spring Energized 


Seal, Enerrin 


Seal, Enerrin 


Threaded Boss 


Threaded Boss 


Shim Stock 


Separator, Gas/Liquid 


Heat Exchanger 


Ferrule, Nut 


Adapter 


Adapter 


Adapter 


Nut, Tubing 


Union, Fitting 


Seamless Tubing 


Regulator, Hydrogen 


Orifice 


Cross 


Adapter 


Adapter 


Tape, PTFE 


e 


Bolt 


Washers, Lock 


Nut 


Insulation 


Brg. Pad Actuator 


Brg. Pad Adjustment 



Alternate Configuration 
1/4-28 X IT 


1/4" to 1/4 FNPT 


1/4" to 1/4 MNPT 


1/4" to 1/2 MNPT 


1/4" 


316 CRES SS 


Cash Acme 


Farmington Eng. 


1 1/4" 


3/16 to 1/4 



Table 1 Recovery System List of Parts 
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Turbocompressor 

AT9801 1 is the turbocompressor assembly drawing. To minimize clearance 
changes due to thermal expansion, it was decided to make the entire 
turbocompressor out of the same material. Additionally, the material had to be 
hydrogen compatible and of sufficient strength to react the rotor centrifugal 
stresses at the initial design point speed of 1.3 million revolutions per minute 
(RPM). Titanium was attractive from a weight, cost, and machinability 
perspective; however, the rotor stresses were marginal. Upon further analysis, it 
was decided to manufacture the turbocompressor from a Titanium alloy (Ti-6AI- 
4V) for the non-rotating parts and Ti-6AI-4V (ELI) STA for the rotor. 

The ELI designates extra low interstitial and the STA is solution treated and 
aged. The solution heat treat and aging increases the material strength to 
provide the required structural integrity while maintaining the coefficient of 
thermal expansion with the non-rotating parts. Note that the detailed finite 
element thermal and structural analysis was performed in Phase I and 
documented in AFAB Technologies, Inc. final report number AT97202-3 dated 
September 15, 1997. 

Accounting for thermal expansion (shrinkage), centrifugal growth, desired and/or 
reasonable operating clearances, and allowance for desired and/or reasonable 
manufacturing tolerances, diametrical fits, radial clearances, and axial 
clearances were calculated and added to the assembly drawing. To eliminate a 
machine-at-assembly operation, the axial distance from the cover flange to the 
compressor shroud was sized so that it would provide the intended operating 
clearance with no shim if the housing was machined to minimum material 
condition (axially) and the remaining parts were machined to maximum material 
condition axially. Since this is an unlikely scenario, the resulting interference will 
be compensated for with the addition of a shim between the cover and main 
housing to achieve the intended axial clearance between the compressor and its 
shroud. 

Each part within the turbocompressor housing is centered and secured radially 
with an interference fit. Additionally, the entire stack is seated axially with the 
cover, which is bolted to the main housing with a four-bolt flange secured with 
%”-28xT’L stainless steel bolts. This joint is sealed with a pressure activated, 
stainless steel spring energized. Teflon jacketed Furon seal (348-129-01-01) 

This seal provides exceptional performance when used with highly polished seal 
glands (typically a 4 - 8 micro-inch surface finish is required for hydrogen 
service). The supply and discharge ports use standard MS9193 bulkhead 
fittings with Enerring seals (4800x-x-2-SPA). 
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The internal parts were designed for assembly into the main housing with either 
an arbor press (a lead-in chamfer is provided in the main housing) or by heating 
and chilling. A fully assembled turbocompressor can be disassembled as 
follows Remove the turbine exhaust fitting and replace with a bolt to constrain 
main housing axially. Remove the compressor inlet fitting and replace with lifting 
eye Remove four flange bolts and pull cover off. Remove compressor end 
journal bearing with jacking screws (each bearing contains three threaded holes 
for jacking screws). Remove rotor. Install lifting eye into threaded thrust bearing 
sleeve. Pull thrust bearing out. Remove turbine end bearing with jacking 
screws. Remove housing constraining bolt and press nozzle ring out of main 
housing with arbor press entering from turbine exhaust. 

The structural and thermal analysis of the turbocompressor assembly was 
completed in Phase I. 

Rotor 

The rotor (AT98003) was designed in four phases. First, the compressor was 
designed to meet cycle performance requirements resulting in identification of 
the compressor geometry, blade coordinates, and rotor speed. Second, the 
journals were sized to match the bearing performance necessary for rotor 
dynamic sta'ril’ty. Third, the turbine was sized to provide the work to drive the 
compressor and overcome bearing and rotor windage losses. Forth, the thrust 
disk was sized to react the axial loads and achieve thrust balance 

The compressor was initially sized using 2-D CFD analysis. The compressor 
operating conditions are summarized in Table 2. Table 3 list the impeller 
geometry. As shown in Table 3, the speed was reduced from 1.3 million RPM to 
1.1 million RPM. Additionally, the initial compressor design has 25 full-length 
blades. The blade coordinates are included on the rotor drawing. The parallel 
wall diffuser and scroll geometry are summarized on Table 4 and Table 5, 
respectively. Figure 1 is a 3-D model of the resulting, initial compressor design. 


Gas 

Hydrogen 

Inlet Pressure 

4.0 atm 

Inlet Temperature 

26 K 

Outlet Pressure 

14.6 atm 

Outlet Temperature 

56.8 K 

Mass Flow Rate 

1 .8 gm/sec 

Gas Power 

0 5968 Kw 


Table 2: Compressor Operating Conditions 
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Direction of Rotation 

Clockwise viewed from inlet 
1,1 million RPM 

^ Speed 

Number of Blades 

25 

Inlet Hub Diameter 

1 .524 mm 

Inlet Shroud Diameter 

2.54 mm 

Inlet Mean Blade Angle 

38 degrees 

Tip Hub/Shroud Diameter 

8.89 mm 

Outlet Blade Angle H 

90 degrees 

Inlet Blade Thickness 

0.0254 mm 

Outlet Blade Height 

0.1016 mm 

Meridional Passage Length 

2.00066 mm 

Blade Passage Length 

2.5146 mm 

Surface Finish 

1 micro-M (40 micro-inches) 

Inlet Shroud Relative Mn 

0.42 

Tip Mn (U2/a0) 

1.1 

Rotor to Stator Clearance 

0.01 mm 


Table 3: Impeller Geometry 


Type 

Vaneless, Parallel Wah 

Width 

,1016 mm 

Inlet Diameter ^ 

9.068 mm 

Outlet Diameter 

13.97 mm 

Surface Finish 

3125 micro-M (12.5 micro-inches) 


Table 4: Diffuser Geometry 


Outlet Diameter 

1.88 

Outlet Mn 

0.17 

Area Variation 

Linear with Angle 

Outlet Diameter of 8 deg. Cone 

2.45 mm 

Gas Velocity 

60 m/sec 

Surface Finish 

0 3125 micro-m (12.5 micro-inches) 


Table 5: Scroll or Collector Geometry 
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Figure 1: Initial Compressor Geometry 


Using this basic design, a range of steady-state operating points were analyzed. 
The pertinent parameters for off-design compressor performance are: polytropic 
efficiency, mass flow, exit temperature, and exit pressure. The tip diameter of 
the compressor establishes the outside diameter of the shaft. 

With the initial compressor design completed, the second phase of the rotor 
design began - journal design. Preliminary bearing designs were completed in 
Phase I. However, the lower design point speed and final compressor tip 
diameter enabled the bearing design to be finalized. A 2-D, bulk-flow CFD 
analysis was performed to identify the stiffness, damping, load capacity, and 
power loss as a function of speed for the minimum and maximum bearing 
clearances. With the bearing design completed, the rotor diameter and length 
have now been established. 

As reported in Phase I, the turbine is a radial inflow, 50 percent reaction type 
developed by Sixsmith at the National Bureau of Standards. Based on the 
compressor and bearing performance described above, the turbine design 
requirements were updated as listed in Table 6. 


Power Required 

800 - 1000 Watt (0.759 - 0.948 BTU/s) 

Rotor Speed 

1.1 Million RPM 

Flowrate 

1.14 gram/sec GH2 (0.0025 Ibm/sec) 

Inlet Pressure 

3 5 atm (51.5 psia) 

Inlet Temperature 

300 K (540 R) 

Expansion Ratio 

3 2 Total-to-Static 

Max, Tip Diameter 

8 -10 mm (0 315-0.394 in) 

Min. Running Clearance 

0 010 mm (0.0004 in) 


Table 6: Turbine Design Requirements 
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A turbine sizing parametric study was conducted based on a total-to-static 
efficiency of 75 percent and optimum tip speed (U/CO) of 0.66 (Note: the turbine 
size varies). A family of curves was generated based on turbine inlet 
temperature and is shown in Figure 2. The resulting efficiency as a function of 
turbine size is shown in Figure 3. 

Based on this analysis, it was determined that the scaled NBS turbine will not 
provide the required performance. The linear scale factor to match tip speed 
produces a flow path with too much annulus area. This is primarily the result in 
the different operating fluid from the NBS turbine (GHe) to the recovery system 
turbine (GH2). The resulting velocities are too low to produce the required high 
specific work. Table 7 compares the major parameters of the NBS baseline 
turbine to the recovery system turbine. Figure 4 shows the change in flow path 
required for the recovery system turbine. The major difference is replacing the 
inlet vane in the NBS turbine with a volute and extending the NBS turbine blades 
with an elliptical profile to improve efficiency. The recovery system turbine will 
then require either supplemental flow rate (as shown in the table) or a higher 
turbine inlet temperature. To retain the self-contained feature of the recovery 
system, the ambient heat exchanger will be coupled with the burn stack to 
provide a higher turbine inlet temperature for the prototype system. 



Figure 2: Turbine Sizing Parametric Study 
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0 0,1 0.2 0.3 0.4 0.5 0.6 0.7 


Isentropic Veiocity Ratio ^ U/CO 


Figure 3: Turbine Efficiency 


Parameter 

NBS Baseline 

Recovery System 

Fluid 

GHe 

GH2 

Power 

0.569 BTU/sec 

0.875 BTU/sec 

RPM 

600,000 

1,100,000 

PT0/PS2 

4.5 

2.73 

U/CO 

0.67 

0.20 

Total-to-Static Efficiency 

0.799 

0.29 

PT-in 

81.44 psia 

138 psia 

TT-in 

54 R 

540 R 

Flowrate 

0.0224 Ibm/sec 

0.0067 Ibm/sec 

ZR 

364.2 BTU/lbm-R 

752.7 BTU/lbm-R 

Or 

1.268 BTU/lbm-R 

3.479 BTU/lbm-R 

Gamma 

1.698 

1 385 

Tip Clearance 

0.0007 in 

0.0007 in 


Table 7; Comparison of NBS Baseline Turbine to Recovery System 
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Figure 4: Comparison of NBS Baseline Turbine Flow Path (top) 
to Recovery System Flow Path (bottom) 
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With the compressor and turbine designs completed, the axial thrust loads were 
calculated to size the thrust disk. Figure 5 and Figure 6 show the pressure 
distribution across the compressor and turbine, respectively. The thrust load for 
the turbine is -11 .9 N and the thrust load for the compressor is 36.76 N for a net 
thrust imbalance of 24.86 N towards the turbine. This is the maximum thrust 
load to be reacted by the thrust disk, thus completing the initial design of the 
rotor 

1 --T - • 



Figure 5: Radial Pressure Distribution Across Compressor 



Kadim ( iTun ) 


Figure 6: Radial Pressure Distribution Across Turbine 


- 13 - 


AfAB [cchnologicj, Inc. 


AT98007-8 


After the initial rotor drawing was completed and manufacturing had begun, the 
manufacturer asked for relief on the compressor geometry Specifically, the 
compressor had a small (insignificant) amount of lean in the blade coordinate 
table and the table identified a tapered blade. The original geometry specified 
an inlet blade thickness of 0.0254 mm (0.001 inches) and a thickness of 0.1016 
mm (0.004 inches) at the exit. To facilitate manufacturing, the geometry was 
modified for a uniform thickness blade of 0.1016 mm (0.004 inches). The impact 
of this increase in blade thickness is a decrease in compressor inlet flow area 
To compensate, the number of blades was reduced from 25 to 24 and the inlet 
hub diameter was decreased from 1 .524 mm to 1 .5 mm. Since it was the intent 
of the compressor design to have no lean in the blade (to minimize bending 
stress) and to have a constant blade thickness, the drawing was revised to 
incorporate these changes with no impact to predicted performance. This 
configuration represents the second impeller design. 

After numerous attempts to have this geometry manufactured, the supplier again 
asked for more relief. The electro-discharge machining (EDM) process could not 
be controlled sufficiently to yield uniform blades. It was felt that a larger 
diameter sinker (plunge) EDM tool, along with thicker blades, would yield an 
acceptable impeller. However, the larger tool and thicker blades also creates an 
interference with the adjacent blade near the inlet where the clearance is only 
0.1016 mm. As a resuH uiis interference would burn the adjacent blade of!' at 
the inlet. 

To improve manufacturability, the impeller was re-designed with thicker blades 
(0.1778 mm) and the 24 blades were re-designed to be 12 full-length blades and 
12 splitters. The use of splitters more than doubled the available tool path at the 
compressor inlet. Additionally, the use of splitters increased the flow path area 
allowing the hub diameter to be increased to the original size of 1.524 mm. 

These changes did not effect the compressor performance at design point; 
however, it did cause a loss in efficiency at some off-design points. This 
configuration represents the third compressor design. 

Again, the supplier expressed concern with the manufacturability of the third 
compressor design. At this point, further re-design was not considered 
necessary as the geometry could be manufactured with either EDM or milling 
The majority of the impeller could be machined with a 0.010 inch end-mill on a 
high-speed spindle (> 20K rpm). Additionally, the inlet could be milled with a 
0.008 inch stub end-mill, provided the inlet blade height could be reduced. 
Therefore, the inlet hub diameter was increased from 1 .524 mm to 1 .778 mm 
and the blade thickness was reduced from 0.1778 mm to 0.175 mm to 
compensate for the reduced flow area on each of the 12 full-length blades and 
12 splitters This represents the forth and final compressor design (and final 
rotor design), which is shown in Figure 7 The performance was determined to 
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be equivalent to the third geometry The structural and thermal analysis of the 
rotor was completed in Phase I for operation at 1 3 million RPM (and over-speed 
capability to 1 45 million RPM); therefore, it was not updated for the lower 
operating speed of 1 1 million RPM 

A rotor dynamic sensitivity study was conducted by Critical Design Analysis, the 
consultant for this project Bearing stiffness, bearing damping, shaft diameter, 
compressor overhung weight, turbine overhung weight, and thrust disk weight 
were all varied to study their effect on the rotor dynamic behavior. The 
consultant’s report is included as Appendix B Additionally, a detailed rotor 
dynamic analysis was conducted on the final rotor design. The rotor dynamic 
stability was predicted to be excellent with the first bending mode occurring well 
above the operating range at 1.5 million rpm and performance relatively 
insensitive to unbalance. The consultant's report is included as Appendix C. 



Figure 7; Final Compressor Geometry 


Journal Bearings 

The initial bearing size analysis was conducted in Phase I. However, this 
analysis was updated based on the latest geometry and operating conditions of 
the rotor, A two-dimensional, bulk-flow CFD analysis was performed to identify 
stiffness, damping, load capacity, and power loss as a function of speed 
throughout the minimum and maximum bearing operating clearance range. The 
resulting bearing geometry was identified on drawing AT98004 This type of 
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bearing was previously designed, analyzed, manufactured, and tested by AFAB 
Technologies. However, the previous bearings were made from Berylco 25 (a 
copper beryllium spring material) while the recovery system bearings are 
Titanium to minimize clearance changes due to thermal expansion. 

The tolerances are relatively loose as the pads have adjustment capability. 
Dowel pins, springs, and set-screws are used to actuate the bearing pads 
radially inward. A 0.010 inch nominally thick web provides radial travel 
capability. A 0.032 inch nominally thick web provides radial stiffness while 
allowing pad rotation for self-stabilization. Three threaded holes are provided 
for jacking screws to facilitate dis-assembly. During startup and shutdown, when 
the rotor speed is insufficient to develop a hydrodynamic film stiffness necessary 
to react the rotor loads, there will be metal-to-metal contact between the 
bearings and rotor. Therefore, several coatings are identified to address this 
potential cyclic life issue. Under normal operating conditions, the rotor is 
supported entirely on a fluid-film. 

AT98010 represents an alternate bearing geometry. The change in materials 
from Berylco 25 to Titanium has resulted in an increase in bearing metal 
stiffness. As a result, there is a concern that the hydrodynamic forces will be 
insufficient to take advantage of the self-stabilization design feature and may 
actually perform as a rigid pad bearir j To address this concern, an alternate 
bearing geometry was designed to move the actuation mechanism from the 
hinge point to the trailing edge of the bearing. This enables adjustment 
capability of the rigid pad bearing, if necessary. The remainder of the bearing 
design is identical. 

Housing 

The main housing has a single counter bore to hold the bearings and volutes, as 
the axial length of the deepest snap is approximately 1 14 inches. The bore 
diameter is the most tightly controlled dimension as it centers all the parts. It 
has a 20 degree lead-in chamfer at the opening to allow parts to be installed with 
an arbor press 

There are three major ports in the main housing: compressor discharge, turbine 
inlet (2), and turbine discharge. The compressor discharge is a single port to 
match the scroll exiting from the parallel wall diffuser. The turbine supply has 
two ports located 180 degrees apart, which feed into the turbine volute supply 
annulus Two ports are used to prevent the formation of a circumferential 
pressure gradient at the inlet of the volute. A large, single port is located on the 
centerline for the turbine exhaust. The through hole diameter of this port was 
increased beyond standard size to provide access to the turbine volute for 
disassembly 
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Additionally, there are three secondary ports; journal bearing adjustment holes 
(6), instrumentation holes (2), and a thrust bearing supply line. Three holes are 
located at each journal bearing to provide access to the pad actuation screws. 
Two instrumentation holes are included (-02 version only) for a proximity probe 
and pressure and temperature sensors. Due to the thin wall necessary to 
provide thermal isolation between the turbine and compressor, there is 
insufficient wall thickness to allow a threaded boss for the thrust bearing supply. 
Therefore, an 0.080 deep spot face was added to permit welding of a 
diameter titanium tube. 

Cover 

The cover, AT98006, is two parts mated into an inseparable assembly. It 
includes the compressor inlet port, parallel wall diffuser, volute, and compressor 
shroud. The inlet port is extended axially so that a length in excess of 4 
diameters (L>4d) is available upstream of the compressor. The parallel wall 
diffuser clearance is set when the volute is mated to the cover. The volute has a 
square cross-section instead of the round cross-section used in the analysis. 

The square cross-section allows it to be machined with an end-mill and has the 
same hydraulic radius as the round cross-section for no impact on performance. 
Th^ compressor shroud profile matches the compressor blade profile A shim 
installed outboard of the seal groove between the cover and housing sets the 
compressor clearance. 

Turbine Volute 

The turbine nozzle ring, AT98007, is also two parts mated into an inseparable 
assembly. Again, the circular cross-sections have been replaced with square for 
ease of manufacture. The inner and outer diameters of the of the volute are 
tightly controlled to provide the proper turbine clearance and housing 
interference fit, respectively. 

Thrust Bearing 

Initially, a spiral groove thrust bearing (SGTB) was sized to react the net axial 
rotor load discussed in the rotor design section. The spiral groove thrust 
bearing design was optimized based on CFD trade studies conducted by 
Broman (1991 ) resulting in the bearing geometry summarized in Table 8. Its 
performance is summarized in Figure 8 and Figure 9. 
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Numhpr nf 


_ ur 

Ratio^n^r Radius to Outerl^TT^ 

uroove Anole 

— im:.:] 

0.8 I 

--ggj!°-°fGl°gg_Dg&th to Operating 

Inner Radius 

Groove Depth 

11 ~ ] 

22 

4.7 mm j 

0.635 mm 1 


Table 8: Spiral Groove Thrust Bearing Geometry 



Figure 8: 


Load Capacity Of SGTB at Nominal Operating Clearance of 0.28 mm 
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Clearance (mm) 


Figure 9: Load Capacity vs. Operating Clearance at 1.1 million RPM 

However, the si7e of this bearing is near structural limits Therefore, the SGTB 
concept was replaced with a hydrostatic thrust bearing. The hydrostatic thrust 
bearing takes advantage of the high-pressure discharge flow from the 
compressor enabling a smaller area to react the same load. This corresponds to 
a smaller thrust disk on the rotor for lower stresses as well as elimination of the 
spiral groove pattern for reduced manufacturing costs. The hydrostatic thrust 
bearing was sized using a 2-D bulk-flow code. Because of the shallow pockets 
required to avoid pneumatic hammer instability, a single annulus was used in 
place of individual pockets. The final thrust bearing geometry is summarized in 
Table 9. 


Number of Orifice Supplies 

8 

Orifice Size 

0.00825 inches 

Inner Diameter 

0.352 inches 

Inner Annulus Diameter 

0.400 inches 

Diameter of Orifice Location 

0.45 inches 

Outer Annulus Diameter 

0.50 inches 

Annulus Depth 

0.00021 inches 


Table 9: Hydrostatic Thrust Bearing Geometry 

The thrust bearing (AT98008) reflects this geometry with additional design 
features The bearing is supplied by an annulus, which receives flow from the 
single V 4 inch tube welded to the outside of the housing. The annulus distributes 
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the flow into the eight 0 0 in(,;ri feeder lines which intersect the 8 orifice lines. 
There is no break-edge at the inner diameter of the thrust bearing to maximize 
performance. The % inch threads in the support allow dis-assembly of the thrust 
bearing. An instrumented configuration (-02) has a 0.188 diameter radial 
through hole for instrumentation access 

Separator 

Every effort was made to find an off-the-shelf gas/liquid separator. However, an 
extensive search did not locate a suitable separator with the performance and 
operating range necessary for the recovery system The residence time to 
condense the droplets was not sufficient. Therefore, a custom designed 
separator was pursued. Several concepts were evaluated. Measuring 
performance of a gas/liquid separator is difficult without sophisticated 
instrumentation (densitometer, which are notoriously inaccurate) to measure the 
quality (fluid loading of the gas phase) of the upstream and downstream flow. 
Therefore, qualitative testing was pursued to further evaluate the leading 
concept. A full-scale Plexiglas separator model was built and tested (see Figure 
10 ). ^ 



Figure 10: Plexiglas Model Separator 
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The Plexiglas separator allowed the fluid generation of the separator to be 
visually tested, and provided invaluable design validation data. The flow is 
accelerated circumferentially as it enters along the outer wall and is accelerated 
as it enters the inner wall. As the velocity increases, the local pressure and 
temperature decrease and the heavier liquid is forced outward where it 
condenses against the inner wall and collects via gravity feed to the bottom. 

The separator operation is shown schematically in Figure 1 1 . 


^'low Outlet 


Flow Inlet 


Inner Wall 


OutletOrain 



Outer Wall 


Metal Mesh Filter 



Flow Inlet 


Liquid Drain 


Figure 11: Separator Operation 

The results of the Plexiglas model facilitated the design of the stainless-steel 
separator identified on AT98009. The final separator is an all metal (passivated 
stainless steel), welded assembly. Like the hydrostatic thrust bearing supply 
line on the housing, %” tubing ports are welded in place for easy installation with 
compression fittings. 

Heat Exchanger 

The hydrogen-hydrogen heat exchanger was sized for parallel, counter flow 
resulting in a minimum of 791 W transfer. All stainless steel, tube-in-tube, 
counter-flow hydrogen-hydrogen heat exchangers (part number 413) were 
purchased from Exergy, Inc., a specialist in miniature heat exchangers. The 
inner tube is %" with 0.035 inch thick walls and the outer tube is V 2 ” with 0.035 
inch thick walls. The inner tube is rated for 4,500 psia operation and the outer 
tube is rated for 2,000 psia service. The pressure drop in water for the inner 
tube is 50 psid at 1 gpm. The pressure drop in water for the outer tube is 25 
psid at 3 gpm The supplied fittings are pipe thread (NPT); therefore, adapters 
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sealed with Teflon tape are used to convert the ports for mating with %” tube 
nuts The heat exchanger includes two 5/16”-18 mounting holes and bolts. 

J-7 Valve 

Several orifice configurations were evaluated for use a Joule-Thompson Valve 
including Bendix bleed orifices, sharp-edged orifices, and converging-diverging 
nozzles. Cv expansion restrictors from Farmington Engineering were selected 
for use as the J-T valve. They are part number CRMS00630xxA, where the xx 
represents the calibrated orifice size in mils (thousandths of an inch). These are 
calibrated, single piece orifices. The design point hydraulic diameter for the J-T 
valve is 0.030 inches. Therefore, six orifice sizes were procured and tested, 
which ranged from 0.025 inches to 0.035 inches. The orifice is installed in a 
standard %" tube fitting (union) by drilling a 0.156 inch counter bore in one side 
of the union 0.170 inches deep. 

Tubing 

The recovery system components are connected with stainless steel seamless 
tubing. AISI 316 stainless steel YT tubing with 0.035 inch thick walls is used 
throughout. All tubing is field routed and cut to fit. All connections are made 
with industry standard compression fittings. 

Fittings 

A variety of industry standard (either Parker or SwageLok) adapters, unions, 
tees, crosses, and fittings are used to connect the various components. 

Frame 

The recovery system is mounted within a metal frame formed by 1 14 inch slotted 
angle and secured with 14 - 28 bolts, lock washers, and nuts. The rectangle 
formed by the outer perimeter of the angle is filled with insulation. The turbine 
supply and exhaust ports situated to protrude outside of the insulated area. 

Final Analysis 

A flow model was created to conduct the steady-state and transient performance 
analysis and is shown below in Figure 12. The results from the steady-state, 
design-point analysis are summarized in Table 10. 
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Figure 12: Liquid Hydrogen Bo.. Off Recovery System Flow Model 


BOUNDARY NODES 


NODE 

P(PSI) 

T(F) 

RHO 

AREAN(IN'' 



3 

214 

-357.4 

0425 

0 

11 

58.8 

-413 

3.9337 

0 

19 

15 

-100 

0.0078 

0 

SOLUTION 





INTERNAL 

NODES 




NODE 

P(PSI) 

PT(PSI) 

T(F) 

Z 1 


4 

213 9944 

213.9976 

5 

213.9842 

213 9874 

6 

213.9724 

213.9755 

7 

213 4991 

213 7808 

8 

212.4321 

212 7145 

9 

211.7316 

211 9171 

0 

59.0023 

86,4716 

2 

54,5441 

55 7241 

3 

46.8711 

48 1499 

4 

34.2414 

36 346 

5 

17 2116 

20 0495 

6 

16,7329 

16 8127 

7 

16 2406 

23 6461 

8 

15 6321 

15 7335 


357 3994 

357 3999 

358 4516 
357 4206 
398 3398 
-398 469 
412 7254 
413.4317 
334 4433 
336.5278 
-26 9939 

-27 024 
-69 4122 
-69 4634 


IQUALITY 


0.9256 
0.9256 
0.9228 
0.9258 
0.3812 
0.3697 
0 5452 
0 5544 

0 9926 
0.9943 
1,0015 
1,0015 

1 0015 
1 0015 


0 425 
0.425 
0 4307 
0 4241 
1,7091 
1.7605 
0.4338 
0 4004 
0 0709 
0.0526 
0 0075 
0.0073 
0.0078 
0 0075 


0.687 
0 6905 
0 
0 
0 
0 
0 
0 
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BRANCH 

K (LBF-SEC^2/ 
(LBM-FT)''2) 

DP 

(PSI) 

FLOWRATE 

(LBM/SEC) 

VELOCITY 

(FT/SEC) 

Re 

M 

107 

4 60E+06 

0.501 

0 004 

78 4 

2.06E+05 

3 59E-02 

103 

4 70E+06 

5 59E-03 

0 0004 

8 28 

2 18E+04 

3 79E-03 

104 

3.20E+06 

1 02E-02 

0 0004 

8 28 

2 18E+04 

3 79E-03 

105 

4 70E+06 

1 19E-02 

0.0004 

8.28 

2 18E+04 

3 79E-03 

106 

1 50E+11 

198 

0.0004 

397 

1 53E-^05 

0.182 

108 

4 60E+06 

1.07 

0.004 

78.6 

2.06E+05 

3.59E-02 

109 

3.00E+06 

0.701 

0 004 

31.7 

2.38E+05 

5 50E-03 

110 

9,00E'^08 

153 

0.004 

380 

8.19E+05 

6 24E-02 

111 

6.30E+06 

0.202 

0.0012 

37.9 

6.34E+04 

2 26E-02 

112 

4.10E+07 

4.46 

0.0027 

159 

1.92E+05 

9 45E-02 

113 

1.00E+08 

7.67 

0.0027 

172 

1.90E+05 

0.105 

114 

1.60E+08 

12.6 

0.0027 

524 

1.64E+05 

0.234 

115 

2.20E+08 

17 

0.0027 

707 

1.67E+05 

0 319 

116 

6 20E+06 

0.479 

0.0027 

315 

1.78E+04 

8.25E-02 

117 

6.40E+06 

0.492 

0.0027 

324 

1.78E+04 

8.49E-02 

118 

6.00E+06 

0.609 

0.0031 

347 

2.20E+04 

9 61E-02 

119 

6.20E+06 

0.632 

0.0031 

361 

2.20E+04 

9.99E-02 


Table 10: Steady-State, Design-Point Results from Flow Model 

A sensitivity analysis was performed for the cycle. The sensitive variable in the 
cycle IS the temperature entering the JT valve. This temperature directly 
controls the amount of liquid produced by the system. Figure 13 shows the 
results of the temperature sensitivity study. Please note that the y-axis scale is 
percent of goal, which is 35 percent. Therefore, 50 percent of goal represents 
50 percent of 35 percent, or 17.5 percent total liquid yield. As you can see from 
the figure, one degree (from 34 to 35 K), represents a drop in liquid production 
from 100 percent of goal to 38 percent. Another increase of one degree (from 
35 to 36 K) further reduces liquid production from 38 percent of goal to 10 
percent. 



4 35 35 5 36 56 5 

Figure 13: Sensitivity of Liquid Production to JT Valve Inlet Temperature 
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The graph underscores that the heat exchanger effectiveness will ultimately 
control the amount of liquid produced by the system. Therefore, the cycle was 
re-evaluated to determine how the desired J-T inlet temperature can be 
achieved more reliably, thus mitigating the temperature sensitivity risk and 
making the overall system more robust. 

Several alternate configurations were considered. The most robust system 
utilizes a second heat exchanger cooled with bypass boil-off. The inclusion of a 
second heat exchanger increases the robustness of the recovery system by 
diverting gaseous hydrogen flow at 4.0 atmospheres and 26 Kelvin directly from 
the boil-off and into the heat exchanger cooling tube, thus bypassing the 
compressor inlet. This diverted, or bypass, flow is then expanded isenthalpically 
through an orifice to approximately 1.1 atmosphers and 22 Kelvin. The bypass 
flow then enters the heat exchanger to aid in cooling the gaseous hydrogen 
exiting the compressor. This configuration is shown schematically in Figure 14 
and the revised flow model is shown in Figure 15. 


from separator 


to j-T valve L 


to turbine hx 


j from compressor exit 


to atmosphere 


./ 



^ 

to compressor inlet 


Figure 14: Schematic of Bypass Flow Configuration 



Figure 15: Flow Model of Bypass Flow Configuration 


An off-design analysis was conducted. The pertinent parameters for off-design 
compressor performance are: polytropic efficiency, mass flow, exit temperature, 
and exit pressure. These parameters were evaluated and are presented in 
Figures 16 through 19 . 



Mass Flowrate (gm/sec) 


Figure 16: Stage Performance of Baseline Hydrogen Compressor 
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Massflow (gm/sec) 


Figure 17: Flow Rate vs. Pressure Ratio of Baseline Conv^ressor 



Massflow (g/scc) 


Figure 18: Compressor Exit Temperature as a Function of Mass Flow 
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Figure 19: Compressor Exit Pressure as a Function of Mass Flow 

The flow model was run for several transient and quasi steady-state points to 
simulate start-up, shut-down, and over-speed conditions. The critical period for 
the turbocompressor is during startup. A successful start involves accelerating 
the shaft as quickly as possible to a speed at which the thrust bearing and tilting 
-pad bearings capacities exet-ed the operating loads. The impulsive startup 
force will be provided by the pressure ratio across the turbine. As the 
temperature feeding the turbine rises, at a fixed pressure ratio, the output power 
will increase and shaft will turn faster. As speed increases, the power 
requirements of the compressor, bearings, and windage increase with speed 
squared. As a result, the system is self-regulating for operating speed. 

Hardware Procurement 


Upon initial completion of the drawings, drawing packages were sent to three 
precision machine shops for quotes. A purchase order was placed with Dynamic 
Engineering, Inc., as they were the low bidder. With the exception of the 
compressor blade geometry, very few manufacturing issues were identified. The 
few issues that were identified involved drawing interpretation or clarification and 
the addition of missing dimensions. 

With the exception of the rotor, most hardware was delivered on schedule or 
within a few months of the promised delivery date, although some hardware was 
delivered incomplete or not at all in accordance with the revised purchase order. 
The numerous re-designs of the compressor to facilitate manufacture delayed 
hardware availability by over a year The design cycle for each of the four 
compressor designs took approximately three months including design, analysis, 
drawing changes, and machining attempts. 
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After the forth and final compressor design, DEI was still unable to machine the 
rotor or find a subcontractor willing to machine it. Despite the inability of a 
precision machine shop with an excellent reputation to produce the part, AFAB 
took receipt of the hardware “as-is" and began investigating alternate suppliers 

The final rotor design was then sent to ten precision machine shops. Two shops 
with extensive turbomachinery experience quoted the job, NREC and Turbocam 
NREC, partnered with Mikino, was selected as the most capable supplier for this 
part size (and they were also the lower bidder). NREC completed the rotors 
within a couple months of schedule, but still beyond the originally scheduled 
completion date for the 24-month contract. Therefore, a 4-month no-cost 
contract extension was requested and approved. The only manufacturing issue 
with the rotor was the recessed radius on each side of the thrust-bearing disk 
was replaced with a conventional radius. 

The Phase II proposal estimated that low volume production costs of the 
turbocompressor would be around $10,000 with a total liquid hydrogen boil-off 
recovery system unit cost of around $40,000. The prototype turbocompressors 
were delivered at approximately $12,000 each with low volume production costs 
estimated at just less than $1 0,000. However, the prototype systems were 
delivered at approximately $20,000 each (hardware and assembly). The final 
production configuration system is now estimated at less than $30,000 each 
(fully assembled and tested), approximately 25 percent less than initially 
estimated. 

Component Tests 

Component testing includes both sub-components and components. The sub- 
components are the bearings and the Joule-Thompson valve. The components 
are comprised of the turbocompressor, heat exchanger, and separator. A data 
acquisition system consisting of 12 thermocouples, 4 pressure transducers and 
a high frequency, AC coupled proximity probe was assembled. The system runs 
on a Pentium II 350 running Windows and LabVIEW software. Temperature 
measurements are resolved to 0.1 degree F and pressure measurements to 0.05 
psi. The AC coupled high frequency data is recorded at 400 KHz. A summary of 
the data acquisition is provided in Table 1 1 below. All instrumentation is fully 
calibrated and traceable to the National Institute for Standards and Testing 
(NIST) 
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1 Manufacturer 

Description 

Omega 

Type “E" (Chromal Constantine) Thermocouples 

Kaman 

KD 2300-0. 5SU Proximity Probe System 

Validyne 

DP-10 Welded SS Transducer 

Validyne 

P2 Transducer 

Validyne 

CD280 Carrier Demodulator (2 Channel) 

National Instruments 

PCI-MIO-16EX I/O Board 

National Instruments 

PCI-MIO-16E1 I/O Board 

i National Instruments 

LabVIEW Software 


Table 11: Data Acquisition System 


Stiffness tests were conducted on the bearings as a function of web thickness. 
The web thickness was rneasured on the minimum and maximum thickness 
bearings and several thickness’ distributed throughout this range. Load and 
deflection measurements were made to calculate the sprint constant as a 
function of web thickness. 

Five orifice sizes, which will be used as the Joule-Thompson valve, were flow 
tested at design point and off-design in air, fioi um, at id hydrogen. Pressure and 
temperature sensors were placed upstream and downstream of the orifice along 
with a flow meter. The data was used to calculate the flow parameter and JT 
factor of the orifice. 

Several turbocompressor builds were assembled and tested. Rotating tests of 
the turbocompressors were attempted in air and helium. Compressor inlet flow 
rate, pressure, and temperature were measured along with discharge pressure 
and temperature, Turbine inlet pressure and temperature, along with flow rate, 
were also measured. The data was to be used to calculate the compressor and 
turbine efficiencies 

Heat exchanger flow tests were conducted at design point and off-design in air, 
helium, and hydrogen. Inlet and discharge pressure and temperature, as well as 
flow rate, were measured for each of the parallel, counter-flow streams. The 
data was used to calculate the heat exchanger efficiency at design and off- 
design conditions 

Finally, qualitative tests were conducted on a Plexiglas model separator and flow 
tests were conducted on a stainless steel separator in air, helium, and hydrogen. 
The separator was tested at design and off-design conditions. 
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System Tests 

Based on the results of the component testing, the optimum recovery system 
components were assembled into a liquid hydrogen boil-off recovery system and 
tests were conducted in air, helium and hydrogen (less turbocompressor) Tests 
were conducted at design point and off-design operating conditions 

Prototype Tests 

The final configuration resulting from the system tests was assembled into a 
prototype recovery system without instrumentation. This system was tested in 
helium prior to delivery to Stennis Space Center. 


- 31 - 


AfAB icchnoloqici, Inc. 


AT98007-8 


Results Obtained 

Provided below are the results obtained for the work performed identified in the 
previous section 

Detail Design 

To validate the integrity of the compressor design, a 3-D heat transfer analysis 
was conducted. The results of this analysis are presented in Figure 20. The 
change in performance from the baseline compressor to the final compressor 
(with splitters) is shown in Figures 21 and 22. 



Temperature 

-216.53 
-220.94 
-225.35 
-229.76 
“234J7 
-238.57 
-242.98 
-247-39 


Figure 20: Thermal Contour Plot (degrees F) of Compressor Blade 
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Figure 21: Stage Performance of Hydrogen Compressor 



Figure 22: Flow Rate vs. Pressure Performance 


The final turbine design resulted in the following recommendations for the 
production turbine: 
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• Vane Ring 

• Provide highly-swirled discharge flow at near sonic speed 

• Exhibit low trailing edge (TE) thickness (<10 percent TE blockage) 

• Have smooth surface finish 

• Radial-Inflow Turbine Wheel 

• Radial star zero-cambered “inducer” blades 

• Axial-turbine-like cambered exducer to provide increased reaction 

• Low trailing edge thickness (<10 percent TE blockage) 

• Have smooth surface finish 

• Tail cone to minimize dump loss and maximize efficiency 

• Exhaust Diffuser 

• Produce low rotor discharge pressure for maximum expansion ratio 

• Maximum practical length and an equivalent cone angle of 2-theta 
(6 degrees) 

These results are presented graphically in Figure 23. Note that the star and 
exducer can be manufactured as separate pieces and later mechanically 
connected or bonded together. The operating conditions and performance for 
the production turbine are shown in Table 12. 



▲ 


0.125” R 

Figure 23: Production Turbine Geometry 
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Fluid 

GH2 

Power 

0.759 BTU/sec(800 W) 

RPM 

1,100,000 

PT0/Ps2 

2 90 

U/CO 

048 

Total-to-Static Efficiency 

0.63 

PT-in 

51.45 psia (3.50 atm) 1 

TT-in 

540 R (300 K) 

Flow Rate 

0.0025 Ibm/sec (0.38 cmm) , 

ZR 

752.7 BTU/lbm-R 

Cp 

3.479 BTU/lbm-R 

Gamma 

1.385 

Tip Clearance 

0.0004 inches 


Table 12: Production Turbine Performance and Operating Conditions 

Performance of the journal bearings was re-calculated based on the final design 
and operating conditions. The resulting performance is shown in Figure 24 
through Figure 28. 



Figure 24: Journal Bearing Stiffness vs. Speed and Clearance 
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Figure 25: Journal Bearing Damping vs. Speed and Clearance 



Figure 26: Journal Bearing Load Capacity vs. Speed and Clearance 



Figure 27: Power Loss vs. Speed and Clearance (Compressor-End Bearing) 
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Figure 28: Power Loss vs. Speed and Clearance (Turbine-End Journal Bearing) 
Final Analysis 

As a result of sensitivity cycle analysis performed, the final cycle includes the 
use of a bypass heat exchanger to increase the robustness of the system. As 
the bypass flow increases, the flow through the compressor decreases. 
Therefore, the compressor efficiency and pressure ratio change according to 
Figure 29 and Figure 30. 


Compressor Efficiency vs. Bypass Flowrate 



0 01 02 03 04 05 06 07 08 09 1 

Bypass Flow(g/s) 

Figure 29: Compressor Efficiency vs. Bypass Flow 
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Compressor Pressure Ratio vs. Bypass Flowrate 



Figure 30: Compressor Pressure Ratio vs. Bypass Flow 


As the efficiency of the compressor decreases, the compressor exit temperature 
increase': by the relation: 

T — X V 

I exit ““ * inlet X r r\ 

where: 


Texit is the compressor exit temperature, 
Tiniet is the compressor inlet temperature, 
PR is the compressor pressure ratio, 

Y is the ratio of specific heats, and 
ec is the compressor efficiency. 


The percent yield of the recovery system is equal to the mass flow of the liquid 
recovered divided by the mass flow of gaseous hydrogen leaving the tank 

Y :: (mcompressor ~ mseparator gas)/nitank “ (1 ~ (hjT~hf)/(hg-hf)) X (mcompressor/rTltank) 

where: 


Y 

is the 

nricompressor 

is the 

rHseparator gas 

is the 

nritank 

is the 

hjT 

is the 

hf 

is the 


is the 


percent liquid yield, 

mass flow through the compressor, 

mass flow of gas leaving the separator, 

mass flow leaving the tank, 

enthalpy of the gas entering the JT valve, 

enthalpy of liquid hydrogen at 4 atmospheres, and 

enthalpy of gaseous hydrogen at 4 atmospheres 
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The enthalpy of the gaseous hydrogen entering the Joule-Thompson valve 
depends upon the enthalpy change possible through the heat exchanger and is 
governed by the relation: 


bypass (^b2 b 1 ) ’^Tomprcssor 


hji = 


h I 

comp ~ *^\'omprcs.sor ^ , ^1 ^ ^ s2 ) 


, h^l ^s"" 
m ( ■ 

‘"compressor ' , , 

hif-h, 


where: 


nibypass is the mass flow diverted from the compressor inlet, 
hbi is the enthalpy of the bypass flow as it enters the heat exchanger, 

hb2 is the enthalpy of the bypass flow as it exits the heat exchanger, 

hsi is the enthalpy of the gas leaving the separator as it enters the 

heat exchanger, and 

hs2 is the enthalpy of the gas leaving the separator as it exits the heat 
exchanger. 

The enthalpies of the cooling gas streams exiting the heat exchanger depend 
upon the effectiveness of the heat exchanger and the maximum heat transfer 
possible through the heat exchanger. The maximum he'=’^ transfer possible for 
each stream of the heat exchanger would occur if the temperature of the cooling 
stream could be raised to the temperature of the hot stream entering the heat 
exchanger. Therefore, the enthalpies of the cooling gas streams exiting the heat 
exchanger are given by: 


bb2 (^”E) X hb1 "T E X hb2max 


and 

hs2 = ( 1 -E) X hsi + Ex hs2max 

where: 

E is the overall heat exchanger effectiveness, 
hb2max is the enthalpy of the bypass cooling flow exiting the heat 
exchanger raised to the temperature of the gas leaving the 
compressor, and 

hs2max is the enthalpy of the separator cooling flow exiting the 
heat exchanger raised to the temperature of the gas 
leaving the compressor. 

Figure 31 plots the relationship between the amount of mass flow diverted from 
the compressor to the overall liquid yield of the system. 


- 39 - 



AfAB fecfmologici, Inc. 


AT98007-8 



Figure 31: Liquid Yield vs. Bypass Flow 

As shown in Figure 31 , the ratio of liquid yield mass flow to boil-off mass flow 
increases with bypass heat exchanger flow rate until the gas approaches vapor 
saturation temperature at the heat exchanger exit. Beyond this point (0.00198 
Ibm/sec), liquid yield decreases. 

In order to calculate the lengths and diameters of the heat exchanger tubing, 
each heat exchanger was divided into five sections. This sectioning was 
necessary because of the significant change in Cp of the hydrogen as it is 
heated or cooled at such low temperatures (Figure 32 and Figure 33). 
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The results for the length of the heat exchangers using 1/8" inner tubing for the 
high pressure flow and 3/16" outer tubing for the low pressure flow are plotted in 
Figure 34 



Figure 34: Heat Exchanger Length vs. Bypass Flow 


Series 1 in the plot corresponds to the length of the bypass flow heat exchanger, 
while series 2 corresponds to the length of the separator gas flow heat 
exchanger. The length of the bypass exchanger varies from 6 to 8 feet and the 
length of the separator gas exchanger varies from approximately 14 to 18 feet 
The OTS heat exchanger is 20 ft long and uses a larger tube, than the analysis 
indicates is necessary. Both of these factors should work to improve 
performance 

Hardware Procurement 


The results of the hardware procurement are the acquired parts, which are 
shown in the following figures. In general, all parts were received to print. 
However, several non-conformances were noted on the assembly log The 
assembly layout shows a sufficiently large diameter at the turbine exhaust port 
through hole in the main housing to facilitate removal of the turbine nozzle ring 
(volute). However, the detail print did not specify the correct through hole 
diameter (the through hole in the main housing was smaller than the through 
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hole in the nozzle ring) thus pir^venting removal of the nozzle ring for 
disassembly Unfortunately tfiis was not discovered until after the first build 
The remaining housings were re-worked at AFAB to the intended dimension 
Additionally the rotor did not have a recessed radius on each side of the thrust 
disk Instead, it had a conventional radius To prevent an interference, a break 
edge was added to the thrust bearing 
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Figure 38: Turbine Nozzle Ring (Volute) 



Figure 39: Cover 
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Component Tests 

The PC-Based, LabVIEW data acquisition system was setup to interface with the 
input/output cards A test control panel was created to display the acquired data 
and activate data logging A sample screen is shown in Figure 42. All the 
instrumentation came with calibration sheets providing traceability to the 
National Institute for Standards and Testing (NIST) Additionally, the 
calibrations were verified by AFAB Technologies to ensure proper operation with 
the data acquisition system A sample calibration curve for the pressure 
transducers is shown in Figure 43, 



Figure 42: Sample Data Acquisition Screen 
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Figure 43: Sample Pressure Transducer Calibration Curve 

The journal bearinn spring rates were measured as a function of web thickness 
The highly non-linear spring rate is theorized to be the result of the complex web 
geometry and possible influence from the EDM re-cast layer formed during 
manufacturing (Figure 44). 



Figure 44: Journal Bearing Web Springrate 
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The orifices were tested for use as both the bypass flow control and the J-T 
valve Various office sizes were tested to ensure a bypass flowrate of 0 001 98 
Ibm/sec at a pressure ratio of 4. The results are shown in Figure 45 through 
Figure 49. This testing determined that a 0.035 inch diameter orifice provides 
the desired flow rate required to maximize liquid yield. 



Figure 45: 0.026 inch Diameter Orifice 
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The JT orifice was tested at varying flow rates, pressures, and temperatures to 
evaluate its performance The results are plotted in Figure 50 


JT Factor = 


hdownstream/h 


upstream 


Where: 


hdownstream i s the enthalpy of the gas downstream of the JT valve, and 
and hupsifeam is the enthalpy of the gas upstream of the JT valve. 



Figure 50: J-T Factor vs. Flow Parameter for 0.035 inch Orifice 

The heat exchangers were tested at different temperatures, pressures, and flow 
rates to determine the effectiveness of each. The results are shown in Figure 
51. The resulting effectiveness for the expected conditions at the compressor 
discharge is 0.9975. 
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Figure 51: Heat Exchanger Efficiency vs. Flow Parameter 
System Tests 

Several system configurations were built and tested. The major differences from 
build-to-build were the turbocompressor fits and clearances, and the J-T and 
bypass orifice sizes. Numerous issues were encountered during 
turbocompressor assembly. 

The turbine nozzle ring (volute) could not be successfully installed in the 
housing with an arbor press. While the part was being pressed in, it would 
deflect in an oil can manner causing the inboard half of the volute to fall to the 
bottom of the housing. Although it could be re-seated, there was no way to 
ensure the circumferential orientation was correct. Therefore, the part was 
installed after it was chilled in liquid nitrogen and the housing was heated with a 
heat gun. By heating and chilling, the part installed properly without issue. 

As each part was installed, axial dimensions were measured to ensure it was 
completely seated prior to installation of the subsequent part. This process 
continued with the addition of the turbine-end journal bearing, thrust bearing, 
and rotor. Once the rotor was installed, a critical dimension was taken from the 
back of the thrust disk to the top of the housing to calculate the shim necessary 
to set the compressor shroud clearance (see Figure 52). After this process, the 
compressor-end journal bearing was installed, followed by the cover. 


- 53 - 



AT98007-8 


AiAH iechnoicqio'. iti'. 



Figure 52: Turbocompressor Rotor Stack 

The initial build was tested in air and the rotor locked-up immediately. The 
turbocompressor was then disassembled where two additional observations 
were noted First, the disassembly thread on the thrust bearing needs to be 
changed to a standard size (from 0 75 inches - 28 to 0.75 inches - 16) 

Second, the disassembly threads on the journal bearing (for jacking screws) 
need to be changed to a standard size (from #8-36 to #8-32) Although the 
ultra-fine threads are desirable for disassembly, they are difficult to acquire. 

Witness marks on the rotor indicated it was bound radially at the turbine nozzle 
ring A review of the design showed a potential radial interference due to 
insufficient feature control of the relative concentricities for the mating parts 
The tolerance stack-up for concentricity provided the opportunity to more than 
exceed the desired build and operating clearances To address this issue with 
the hardware we had each part was inspected and selected based on its fit with 
the mating part Virtually every part selected was honed and polished to 
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Optimize its fit and clearance to the extent possible The remaining 
turbocompressor builds used this process for hardware selection and were built 
to the previously described build process Unfortunately, the polishing could not 
overcome the concentricity tolerance issue Each test resulted in immediate 
rotor seizure (although the remaining rotors seized in the journal bearings rather 
than the turbine volute). Turbocompressors were tested both by themselves and 
as part of the system tests A system test setup is shown in Figure 53. The 
yellow rectangles are Type "E” (chromal Constantine) thermocouples and the 
cylinders are pressure transducers. 



Figure 53: System Test Setup 


Figure 54 shows the effect of heat exchanger effectiveness on the liquid yield of 
the system. The measured heat exchanger effectiveness results in a liquid yield 
of slightly more than 31 .2 percent. However, because the JT factor of the orifice 
is less than one, the liquid yield of the system is lowered to 31 percent as shown 
in Figure 55 
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Prototype Tests 

The test results from the system testing were used to finalize the orifice sizes for 
the bypass flow control and Joule-Thompson valve A final build of the 
turbocompressor was attempted and the prototype system (same as the system 
test configuration without instrumentation) was assembled. The prototype was 
tested in helium and again, the rotor locked-up immediately. Therefore, the 
prototype was not insulated, thus allowing immediate access to the 
turbocompressor for subsequent re-work The prototype, along with the 
remaining hardware, was packaged and sent to Stennis Space Center (see 
Figure 56). 



Figure 56: Prototype System 
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Technical Merit and Feasibility Assessment 

The liquid hydrogen boil-off recovery system has significant technical merit It is 
based on a modified Linde-Hampson cycle used previously to liquefy gases 
Sub-component, component, and system level tests indicate 31 percent of the 
boil-off from a liquid hydrogen storage tank can be re-liquified. 

Two technical challenges were discovered as a result of the program: 
manufacturability of the compressor blades, and a rotor-to-stator clearance 
problem in the turbocompressor as a result of a concentricity tolerance stack-up 
The compressor blade producibility challenge was overcome during the project 
and the rotors were successfully fabricated. The rotor-to-stator 
clearance/concentricity issue can be resolved by re-designing the parts to be 
line-bored in the housing, thus eliminating the tolerance stack-up. 

However, if the technical challenges with the turbocompressor cannot be 
overcome, the turbocompressor can be replaced with an alternate means of 
compression. When it appeared the rotors could not be manufactured due to 
the compressor blade geometry, AFAB Technologies identified an electrically 
driven, cryogenic piston compressor which could be substituted into the system: 
Fluitron model AI-15-NL. This pump is roughly equivalent in cost and 
oerformance of the turbocompressor and thus will not impact the performance of 
me liquid hydrogen boil-off recovery system. 

In low volume production, the system can be fabricated, assembled, tested, and 
delivered for less than $30K per unit. A single 22,000 gallon low pressure liquid 
hydrogen storage tank loses $20,000 per month due to boil-off. The liquid 
hydrogen boil-off recovery system would save $7,000 per month based on a 
volumetric efficiency of 35 percent. Therefore, the system would pay for itself in 
approximately five months of service. 
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Appendix B 

Consultant’s Report - Sensitivity Study 


Rotordynamic Sensitivity Study 

To analyze the effects of design changes or prediction errors in the dynamic 
characteristics of components, a sensitivity study was performed. The following 
parameters were varied to study their effect on the rotordynamic behavior of the 
rotor 

1 ) Bearing stiffness 

2) Bearing damping 

3) Shaft diameter 

4) Compressor overhung weight 

5) Turbine overhung weight 

6) Thrust disk weight 



7) Bearing Span 

To analyze the effects on the rotor behavior of the design parameters, the undamped 
critical speeds are used to demonstrate changes in the frequency of the modes It is 
desirable to maintain a separation margin between the operating speed range and 
the location of the critical speeds. The baseline design from Phase I provided 
separation margins in excess of 50% (20% is deemed marginal.) The response of 
the rotor to unbalance will be used to quantify the effects of changes in the bearing 
coefficients. The results of the study are summarized below. 

Baseline 

Three levels of stiffness were selected corresponding to the approximate magnitudes 
of the predicted bearing stiffness at minimum, maximum and nominal clearance. 
These are shown below along with the first three critical speeds. 


Beating 

Baseline 

Stiffness (Ibf/in) 

Critical Speed (rpm) 

1 

2 

3 

10,000 

207,735 

227,252 

2,285,921 

50,000 

462,069 

505,440 

2,289,703 

250,000 

1,006,321 

1,098,858 

2,312,605 


The undamped critical speed map and mode shapes are graphically shown below for 
a bearing stiffness of 50,000 Ibf/in. 


Liquid Syut^ro 

Pflromelw Study 

June 6. 19^8 CCA 


Rotor 1 Brq ^ 


(J fc- 



1 r e a r S ^ i M ne ss ( b ‘ n ) 













N - 2280703. RPM j 



STIFFNFSS CASF F 



LB/IN 

LH • IN 

K 1 

1 50000. 

OGOO 

K 1 

2 50000. 

.OGOO 

KC 

1 10000000 

.OGOO 

KC 

2 10000000 

.OGOO 


8 

? 

S 


Liquid Hydrogen Boil- off Recovery System 

Pororri^ter 5en3HK/ity Study 

Jur>e 6. 1998 CDA 


The critical speed map plots the first four modes for a range of bearing stiffness from 
1 000 to 250,000 Ibf/in at the two bearing locations. The mode shapes illustrate that 
the first two modes are rigid body modes assodated with the overhung compressor 
and turbine. (The two shafts shown in the mode shape represent the actual rotor 
displacement and a shorter dummy rotor required by the analytical method.) The 
third mode is the first bending mode of the rotor. The bending modes are highly 
sensitive to unbalance in the rotor. These produce large bearing reactions and shaft 
vibrations. These modes require separation from the operating range to ensure 
smooth trouble free operation. 

The frequencies for the critical speeds will be used as the baseline numbers and are 
based on the Phase I layout of the rotor. While used as a baseline, the critical speed 
map does show an important consequence of the current rotor design. The bending 
mode is located far above the operating speed and is insensitive to the bearing 
stiffness. At greater stiffness levels, the two rigid body modes approach the 
operating speed of 1 ,300,000 rpm. 



Shaft Diameter 


The shaft diameter was varied by ±25% of its current value. The critical speed 
calculation was carried out for each case and the results are tabulated below as a 
frequency and a percentage change over the baseline. 


Bearing 

Stiffness (Ibf/in) 

Shaft Diameter +25% 
Critical Speed (rpm) 

1 

2 

3 

10,000 

173,485 

189,832 

2,947,682 

50,000 

386,946 

423,381 

2,949,694 

250,000 

854,415 

934,364 

2,960,506 

Bearing 

Stiffness (Ibf/in) 

Shaft Diameter -25% 
Critical Speed (rpm) 

10,000 

252,705 

286,888 

1,573,624 

50,000 

554,211 

6j2,967 

1,585,708 

250,000 

1,109,307 

1,308,506 

1,694,537 


What is apparent from the changes in the shaft diameter is that an increased 
diameter increases the frequency of the bending mode (3"^) and decreases the 
frequency of the rigid body modes (1 & 2'^). Since the bending mode is strongly 
affected by the shaft stiffness, a larger diameter shaft will have bending modes with 
higher frequencies. Since the absolute change in frequency was equal for an 
increase or decrease in diameter, only detrimental changes will be examined for the 
remaining design parameters. 

Bearing Span, Thrust Disk Size, Compressor and Turbine Weight 

The effects of increasing the bearing span and thrust disk size by 25% and the 
compressor and turbine weight by 100% were examined next. The results are 
tabulated below: 















Bearing 

Stiffness (Ibf/in) 

Bearing Soan +25% 
Critical Speed (rpm) 

1 

2 

3 

10,000 

198,213 

225,116 

1,746,143 

50,000 

439,857 

499,760 

1,755,880 

250,000 

942,618 

1,072,611 

1,822,723 

Bearing 

Stiffness (Ibf/in) 

Thmst Disk Size +25% 
Critical Speed (rpm) 

10,000 

205,807 

225,129 


50,000 

457,811 


2,244,782 

250,000 

997,541 

1,081,100 

2,275,876 

Bearing 

Stiffness (luf/in) 

Comoressor Weight +100% 
Critical Speed (rpm) 

10,000 

136,330 

208,623 

1,079,786 

50,000 

287,899 

464,168 

1,134,044 

250,000 

501,851 

1,012,338 

1,393,301 

Bearing 

Stiffness (Ibf/in) 

Turbine Weight +100% 
Critical Speed (rpm) 

10,000 

125,602 

223,984 

1,706,844 

50,000 

273,885 

498,138 

1,740,010 


250,000 


544,636 


1,082,898 


1,706,844 
1,740,010 
1 ,902,393 

















The inaease in bearing span, possibly necessitated by an increase in bearing length, 
thrust disk enlargement or spatial restrictions on the exterior of the case, has a large 
impact on the location of the bending mode. Increased bearing span, i.e. shaft 
length, has an inverse relationship with the shaft stiffness. As with the shaft 
diameter, a decrease in shaft stiffness will lower the frequency of the bending mode 
In this case, a 25% increase in span has dropped the bending mode by an equivalent 
amount. As can be seen in the table, the rigid body modes are not significantly 
affected by the increase in bearing span. 

To examine thrust disk changes due to performance variations in either the 
compressor, turbine, or the thrust disk itself, a 25% increase in the diameter and 
thickness was analyzed. The results indicate that this magnitude of change does not 
affect the dynamic behavior of the rotor. Frequencies of the rigid and bending modes 
changed only slightly. 

Performance variations, due to inaccurate predictions, machining errors, or changes 
in inlet/exit conditions, will also affect the compressor and turbine design. One of the 
impacts of these design changes is an inaease the overhung weight at the shaft 
ends. A study was performed where the component weight at each end of the shaft 
was increased by 100%. While this may seem a large increase, changes to the 
overhung length in shaft and, not only, design changes to the turbine and 
compressor, will add weight in a similar fashion. Additionally, the results can be 
accurately interpolated for fractional changes in weight. 

The increase in compressor weight had the greatest impact on the location of the 
critical speeds. Of great interest is the drop in frequency of the bending mode. For 
the increase in compressor weight noted, the bending mode has fallen below the 
operating speed of 1 ,300,000 rpm. This would impair operation of the rotor system. 

The inaease in turbine weight had a lesser impact on the frequency of the bending 
mode. The table shows a deaease of only 1 0% in frequency. The differences are 
due in part to the smaller overhang of the turbine versus the compressor. 

Bearing Stiffness & Damping 

Since the effects due to variations in the bearing damping can not be measured 
using an undamped analysis, the rotor’s response to unbalance will be used. This is 
a forced response analysis and incorporates the bearing stiffness and damping 
characteristics. The predicted response at the turbine, mid-span and compressor will 
be used as the basis for comparison. The peak response and response at 1 ,300,000 
rpm speed will be monitored for changes. Stiffness will be varied by the predicted 
curves for minimum, nominal and maximum clearance. Damping levels will be 
multiplied by % and 2. A fixed level of unbalance placed at the turbine and 
compressor locations is used throughout the analysis. The magnitude of the 
unbalance is not meant as representative of expected levels, but merely selected to 
produce sufficient response levels for the analysis. 



An additional measure of behavioral change for variations in damping is the stability 
of the rotor/bearing system. The log decrement is used to represent the system 
stability level. A positive log decrement indicates stable system. 



Baseline 

Location 

Peak Response 
(mils) 

Frequency (rpm) 

Response @ 
1,300,000 rpm 

Turbine 

2,34 

320,000 

1.88 

Mid-Span 

0.45 

u 

0.48 

Compressor 

2.61 

a 

0,25 


Stiffness (+) 

Turbine 

7.29 

640,000 

3.86 

Mid-Span 

2.45 

710,000 

1.04 

Compressor 

9.01 

710,000 

2.13 


Stiffness (-) 

Turbine 

1.57 

580,000 

1.55 

Mid-Span 

0.29 

310,000 

0.40 

Compressor 

1.62 

310,000 

0.39 


Damping (1/2x) 

Turbine 

2.91 

230,000 

1.86 

Mid-Span 

0.62 

a 

0.48 

Compressor 

3.43 

a 

0.13 


Turbine 

Damping (2X) 

1.90 

600,000 

1.86 

Mid-Span 

0.36 

470,000 

0.47 

Compressor 

1.93 

a 

0.49 
























Damping Levels 

Log Decrement of Critical Speed 

1 

2 

3 

Baseline 

0.697 

0.693 

0.012 

’/2X 

0.343 

0.337 

0.0051 

2x 

1.43 

1.43 

0.025 


In the undamped analysis, increased stiffness due to the tighter clearance condition 
increased the frequency of the rigid body modes. The unbalance response analysis 
also predicted an inaease in the peak response speed, roughly doubled. In addition, 
the response due to the unbalance also inaeases in some cases 3-4x the baseline 
values. In contrast, the reduced stiffness case corresponding to maximum 
clearance, has a minimal effect on the frequency of the peak response and lowers 
the response levels by 30-40%. These results indicate that at the higher stiffness 
levels, the bearing damping is less effective in controlling the rotor response at the 
turbine and compressor than at the lower stiffness levels. 

In dynamic systems, damping is effective in reducing the vibration levels at speeds 
near or corresponding to critical speeds. Surprisingly, at speeds removed from 
critical speeds, increased damping levels will increase the response level of the rotor. 
This can be seen in the table. At the operating speed of 1 ,300,000 rpm, lower 
response levels are achieved at the compressor end with smaller damping levels as 
expected. However and more importantly, the peak response magnitude drops with 
increasing damping. Since this represents the maximum predicted vibration, it is the 
critical factor in determining the operability of the rotor. 

Using an arbitrary level of destabilizing force in the rotor, stability of the rotor system 
improves as the damping is increased. This is as expected and indicates that 
excessive damping levels have not been reached. (Note: As with stiffness, 
unnecessarily large levels of damping will decrease the effectiveness of the damping 
in controlling the response of the turbine and compressor shaft ends to unbalance 
and improving the rotordynamic stability of the system.) 

Recommendations 

The following factors should be maximized to improve the dynamic behavior of the 
rotor/bearing system: 


1 










• Shaft Diameter - Increases in the shaft diameter lower the frequency of the rigid 
body modes and raise the bending modes. This enlarges the effective separation 
region in which trouble free operation can be expected. 

• Bearing Damping - The reduction of the peak response levels and the increase 
in rotor stability indicate that optimal damping levels have not been reached. 
Bearing damping levels can usually be increased by lengthening the bearing. 
However, this does raise the sensitivity of the bearing stiffness to clearance and 
increase the bearing span, both detrimental effects. 

The following are design items or details that should be minimized during the final 

design stages of the system: 

• Bearing Span - Increases in the bearing span will lower the frequency of the first 
bending mode. This could place restrictions on the maximum safe operating 
speed. 

• Bearing Stiffness - The study indicates that stiffness levels from 10,000 to 50,000 
Ibf/in are more desirable for the present rotor configuration than 250,000 Ibf/in. 

The sensitivity of the stiffness to clearance changes can be reduced by 
shortening the bearing length, however this has a negative impact on the 
damping (see above.) 

• Overhung Weight - The overhung '"~'ght (and moment) at the shaft ends should 
be kept as small as possible. This is espedally true at the compressor end. The 
location of the bending mode is strongly dependent on the overhung weight at 
this location. 

• Thrust Disk - While only a minor effect on the dynamic behavior, drastically 
increasing the size of the thrust disk will add weight to the center of the rotor. 

This will lower the frequency of the bending mode and may eventually encroach 
upon the operating speed. 
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Self-Acting, Liquid Hydrogen 
Boil-Off Recovery Sydem 

Final Turbocompressor Configuration 
Rotordynamic Analysis 

Introduction 

The final configuration of the turbocompressor rotor was provided by AFAB 
Technologies, Inc. The dimensions of the shaft were taken from drawing 
AT98003, 10/1/98. As with previous configurations, the dynamic behavior 
and integrity of the design is measured by performing a rotordynamic 
analysis. The analysis includes an undamped critical speed study, 
unbalance response calculation and a stability prediction. 

The shaft material is titanium. Correspondingly, the density and Young’s 
modulus used were 0.167 Ibm/in^ and 14.8 x 10® Ibf/inl The pitch diameter 
was u -ef; to model the mass and inertia of the compressor and turbine 
blading. Since the blading occurs at the shaft ends, shaft stiffness 
contributions of these components are neglible. Maximum operating speed 
is defined as 1 ,1 00,000 rpm. 

The rotor/bearing system of the turbocompressor is displayed on Figure 1). 
The shaft is divided into 16 mass stations, each having a length-to-diameter 
ratio of less than 1 .0. The beginning of each mass station contains the mass 
and inertia effects the surrounding shaft. Shaft stiffness is modeled in 
massless elements connecting the mass stations. 





Relative Displacement (Dim) 


Grid spacing shown is 0. 1 The overall length of the rotor is roughly 1 .6'' and the 
bearing span is *0.87". Bearing positions are indicated by the vertical spring 
schematics. 


Undamped Critical Speed Analysis 

The undamped analysis is used to give a general idea of the locations of the critical 
speeds and mode shapes of the turbocompressor. In this analysis, the bearing 
support stiffness is assumed and varied from 5000 to 500,000 Ibf/in. Damping is 
neglected. 

The first four mode shapes and corresponding frequencies for a support stiffness of 
5000 Ibf/in are shown on Figure 2). The first two modes are considered rigid body 
modes showing little bending in the shaft. Given this support stiffness, these would 
be expected to occur near 128,000 rpm in the operating range. The bearings are 
expected to absorb most of the potential energy of these modes since the shaft is 
behaving as a rigid body. The first two undamped bending modes occur at speeds 
above 1,500,000 rpm. These modes should be avoided during operation of rotating 
equipment. Since the energy produced by rotating forces is absorbed by the shaft 
bending, the bearing have little control on the amplitudes of vibration. 


— 

I AFAB Technologies, Inc. 
Turbocompressor Rotor Model 
Final Configuration June 11, 1999 
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Figure 2) Undamped Mode Shapes at 5000 Ibf/in 
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The analysis is repeated using a support stiffness of 500,000 Ibf/in. In this 
instance, the bearing stiffness is now a larger percentage of the bending 
stiffness of the shaft. As expected, the first two mode shapes of the shaft 
now show some bending. Since the bearings are assumed stiffer by two 
orders of magnitude, shaft bending is now absorbing some of the potential 
energy Additionally, the increased support stiffness has moved the 
frequency of these modes to =1 ,200,000 rpm. To avoid operating near these 
critical speeds, bearing stiffness should be kept less than 500,000 Ibf/in. The 
first two bending modes are still above 1,500,000 rpm and show little effects 
of the increase in support stiffness. This is due in part to the fact that the 
node points of these modes are occurring near the bearing locations. 
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Undamped Critical Speed Map 



Figure 4) Undamped Critical Speed Map 


Unbalance Response Anaiysj^^^^ 

The unbalance response analysis is used to determine how the rotor will react to 
unbalances caused by manufacturing deviations, machining tolerances, and wear. 
Three unbalance cases are examined in this study; compressor, turbine, and thrust 
disk unbalance. The turbine and compressor unbalance distributions should excite 
the first two modes predicted by the undamped analysis. The thrust disk unbalance 
will excite the first bending mode since it is located at the rotor center near the 
maximum deflection of that mode. Unbalance magnitudes of 0.00001 oz-in are used. 
However, the response is linear and can be directly scaled with unbalance levels 

The predicted bearing behavior versus speed is used in this study. Both stiffness 
and damping coefficients are incorporated in the calculations. As with the undamped 
analysis, the bearing stiffness determines the location of the critical speeds. The 
damping will determine to what extent the rotor will respond to the unbalance 
excitation. (As noted before, poor placement of the critical speeds, either by the shaft 
design or bearing stiffness, may hamper the effectiveness of the bearing damping in 
controlling the motion of the shaft.) The bearing coefficients used are contained in 
Table 1) 



Speed (rpm) 

Stiffness (Ibf/in) 

Damping (Ibf-sec/in) 

500,000 

18,000 

0 18 

1,000.000 

40,000 

0.25 

1,500,000 

70,000 

0.37 


Table 1) Bearing Stiffness and Damping vs. Speed 


Since the bearing coefficients are predicted to be symmetric (X dir = Y dir), the rotor 
response in the X- & Y-directions will be identical as well. The rotor response will be 
presented as the amplitude and phase angle at the compressor, thrust disk and 
turbine locations. Critical speeds are normally identified as peaks in the response 
curves. Additionally, rapid changes in the phase angle are also usually connected 
with the presence of a critical speed. 

The first case presented is the rotor response to turbine unbalance. Figures 5-7) plot 
the response of the three stations. Notice that no peaks are evident in the response 
especially near the operating speed of 1,100,000 rpm. It is likely that the first two 
modes are occurring near 200,000 rpm. However, the bearings are effective in 
damping out the response of these rigid body modes. (The rapid phase change 
evident in the turbine response is due to the transition of the mode shape from 
conical to bending. The turbine response passes through zero at this point.) 
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Figure 5) Turbine Unbalance - Compressor Response 
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Figure 6) Turbine Unbalance - Thrust Disk Response 
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Figure 7) Turbine Unbalance - Turbine Response 



The second case presented is the rotor response to compressor unbalance. Figures 
8-10) plot the response of the three stations. Once again, no peaks are evident in 
the response especially near the operating speed of 1,100,000 rpm As before, it is 
likely that the first two modes are occurring near 200,000 rpm. {The rapid phase 
change is evident in the compressor response as it passes through zero at this 
point.) 

Finally, the response due to thrust disk or mid-span unbalance is presented. This 
distribution is intended to excite the third mode or first bending mode of the 
turbocompressor. It is predicted to occur above 1,500,000 rpm so the peak response 
will not be seen. Figures 11-13) contain the predicted response at the same three 
stations as before. No peaks are evident. However, the rapid rise in the response 
near the 1 ,500,000 rpm speed is indicative of the close proximity to a bending mode 
of the system. (The rapid phase change is evident in the thrust disk response as it 
passes through zero at this point.) 
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Figure 8) Compressor Unbalance - Compressor Response 
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Figure 9) Compressor Unbalance - Thrust Disk Response 


AFAB Technologies, Inc 

Turtxjcom pressor Final Configuration - June 1 1 , 1999 
Unbalance Response Analysis 



URoloi St 
URolor St 
URoloi St 
URotor St 


15 X Mug 

16 X Pn<*-i* 
10 Y-Msg 
16 Y Plv»i* 


! s 



Speed (rpm) 


I 


Unbalance information 

SI 1 Co<npt*%so( UnbatAOi.* 

Q 00001 oz tr> 


Stations Descnptions 

St 1Coftipt«iSo» 

St 6Cotnpf*ssof (3tg 
St 9Thrusl Dtsf 
St 1 3 1 urbn* rtf g 
S< ISIurbvi* 


Figure 10) Compressor Unbalance - Turbine Response 
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Figure II) Tlimst Disk Unbalance - Compressor Response 
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Figure 13) Thrust Disk Unbalance - Turbine Response 


S^bility Anaiysls 

The rotor stability is critical to enable operation at full speed. Rotor instability can be 
caused by many factors including the aerodynamic forces in compressors and 
turbines. Rotor instability normally appears as a re-excitation of the first critical 
speed. This occurs while the rotor is operating at speeds greater than twice the 
frequency of the first mode. The vibration can become unbounded and lead to 
severe rotor damage and large transmitted forces. 

Stability is measured by the logarithmic decrement (log dec). A stable system will 
have vibration levels that will return to normal levels after experiencing a pertubation 
(positive log dec). An unstable system will have levels that grow nearly unbounded 
following the perturbation (negative log dec). 

Figure 14) illustrates the rotor model used in the stability analysis. Two additional 
springs are shown at the compressor and turbine locations. These represent the 
destabilizing forces placed at these locations. The destabilizing force is a cross- 
coupled stiffness term. Cross-coupling causes a movement in a direction 
perpendicular to the applied force. A cross-coupled stiffness term induces whirl or 
the tendency towards unbounded growth in the vibrations. The bearing damping 
counteracts these terms. A magnitude of 300 Ibf/in is placed at each location. This is 
considerably larger than expected for this size of equipment operating in hydrogen. 

The results of the stability analysis are shown on Table 2). The operating speed 
used in the analysis is 1 ,100,000 rpm. Maximum speed represents the worst case 
scenario for the analysis. 
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